A(1.00(}'303 


Ship  Structural  Response  Analysis 
Spectra  and  Statistics 


B.  Phelps 


^:^ELECiEl  ^ 
^  NOV  0  9  1995  I  \ 


f^^isrrM.us'jTiyrs  STJiT£M£ici  k 


j  Approvsd  tea  gtiblie  release 
^  Cfistiibmioo  Ualizaitwl 


Approved  for  Public  Release 


©  Commonwealth  of  Austral!^ 


Department ^  of  Defence 

DEFENCE  SCIENCE  AND  TECHNOLOGY  ORGANISATION 


Ship  Structural  Response  Analysis: 
Spectra  and  Statistics 

B.P.  Phelps 


Ship  Structures  and  Materials  Division 
Aeronautical  and  Maritime  Research  Laboratory 

DSTO-TR-0183 


ABSTRACT 

Wave  loading  of  ships  and  maritime  structures  is  a  random  process  and  the  response 
of  these  structures  is  itself  a  random  process.  Much  work  has  been  undertaken  to 
better  understand  the  nature  of  waves  and  a  number  of  existing  techniques  are 
described.  Statistical  and  spectral  analysis  techniques  may  be  used  to  quantify  the 
wave  loads  and  the  corresponding  structural  responses.  The  relationship  between  the 
input  load  and  response  may  then  be  determined  through  the  use  of  response 
amplitude  operators  and  the  response  of  the  structure  to  predicted  lifetime  extreme 
loads  and  fatigue  loading  then  calculated.  Standard  wave  spectra  and  directional 
spreading  factors  are  used  to  enable  the  seaway  to  be  described  mathematically  and 
several  standard  forms  are  discussed.  The  use  of  statistical  methods  enables  the 
seaway  to  be  described  in  terms  of  a  limited  number  of  parameters  from  which  short 
term  and  long  term  probability  distributions  may  be  obtained.  These  distributions 
may  then  be  used  to  enable  extreme  lifetime  loads  to  be  estimated. 
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Ship  Structural  Response  Analysis: 
Spectra  and  Statistics 


Executive  Summary 

Because  of  the  essentially  random  nature  of  ocean  waves,  and  hence  ship  structural 
responses  to  the  loads  imposed  by  waves,  statistical  and  spectral  analysis  techniques 
are  commonly  used  to  quantify  ocean  waves  and  the  corresponding  structural 
responses.  Details  of  these  methods  and  the  manner  in  which  they  may  be  applied  are 
stunmarised  in  this  report. 

The  report  summarises  many  of  the  methods  used  and  the  way  in  which  they  may  be 
applied  in  practice.  For  example,  the  use  of  wave  power  density  spectrum  to  quantify 
ocean  waves,  the  relationship  between  actual  and  encoimter  wave  frequencies,  the 
transformation  of  the  spectrum  between  the  two  frequency  domains  and  the  use  of  the 
directional  energy  spectrum  to  account  for  wave  spreading  are  described.  For  design 
purposes  a  number  of  idealised  wave  spectrum  formulations  are  available  and  the 
applicability  of  the  various  forms  is  discussed.  The  use  of  Response  Amplitude 
Operators  (RAOs),  which  provide  the  ratio  of  the  response  spectrum  to  the  input 
spectrum,  are  discussed  and  the  need  to  pay  particular  attention  to  the  relationship 
between  encounter  frequency  and  actual  wave  frequency  is  highlighted. 

The  long  term  response  of  a  ship  to  wave  loading  over  a  wide  range  of  operating  and 
environmental  conditions  determines  the  fatigue  performance  of  the  vessel  and  the 
report  describes  the  methods  used  to  make  long  term  fatigue  predictions  based  upon 
short  term  observations.  The  use  of  short  and  long  term  statistical  methods  is 
presented  and  the  applicability  of  particular  methods  is  discussed. 
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1.  Introduction 

The  in-service  loading  of  a  ship  may  be  considered  to  be  broadly  made  up  of  two  parts, 
the  still  water  loading  and  the  wave  induced  loading^.  Still  water  loading  is  a  function 
of  the  difference  between  the  mass  and  buoyancy  distributions  along  the  length  of  the 
ship  and  so  is  influenced  by  the  cargo,  fuel  and  stores  condition  and  although 
essentially  random  in  nature,  for  most  purposes  it  is  normally  considered 
deterministic.  Wave  loading  is  a  truly  random  process  as  it  is  not  possible  to  predict  at 
any  given  time  the  height  of  a  wave  relative  to  the  ship. 

Traditionally  analysis  of  wave  loading  has  been  reduced  to  a  quasi-static  condition 
with  the  ship  balanced  on  the  crest  or  between  two  crests  of  a  "standard  wave" 
representing  an  extreme  loading  condition,  and  the  load  is  then  superimposed  upon 
the  still  water  load.  Although  this  method  has  worked  well  for  many  years  it  ignores 
the  random  natvue  of  the  waves,  their  dynamic  effects  and  the  influence  that  this  has 
on  the  total  loading  and  fatigue  life  of  the  structure.  A  better  understanding  of  a  ship's 
response  to  wave  loading  is  therefore  vital  in  the  overall  assessment  of  hull  structural 
performance.  Wave  loading  causes  a  dynamic  ship  response  which  may  be  considered 
made  up  of  a  whole-ship  rigid-body  response  and  a  series  of  deformation  responses 
involving  bending  and  torsion  of  die  hull  girder.  Structural  dynamic  theories  have 
therefore  been  developed  whereby  the  ship  is  considered  to  respond  as  the 
superposition  of  responses  of  increasingly  higher  order  modes  of  vibration.  While 
such  theories  are  not  the  subject  of  this  report  the  idea  of  the  structmal  resporise  being 
the  superposition  of  responses  at  different  frequencies  is  inherent  in  the  methods 
presented. 

Because  waves  and  the  subsequent  ship's  responses  are  random  in  nature  they  are  able 
to  be  analysed  using  statistical  methods,  and  much  work  has  been  rmdertaken  to  better 
understand  the  nature  of  waves  through  their  spectral  and  statistical  distributions  and 
the  corresponding  responses  of  the  ship.  T5rpically  the  seaway  is  described  by  wave 
height  spectra  and  linear  relationships  between  the  wave  height  and  the  ship's 
response  spectra  may  be  developed.  These  relationships  are  often  called  Response 
Amplitude  Operators  and  their  linear  nature  means  that  their  application  is  restricted 
to  fairly  mild  and  uniform  sea  conditions,  they  do  not  differentiate  between  hogging 
and  sagging  responses  and  effects  such  as  breaking  waves  and  slamming  are  not 
adequately  described.  There  still  remains  therefore  much  work  to  be  done  before  the 
relationship  between  ship  and  waves  is  fully  imderstood. 

The  development  of  a  broadly  based  expertise  in  available  structural  analysis  methods 
is  considered  essential  to  maintaining  ^e  ability  for  the  RAN  to  assess  the  structural 
performance  of  existing  and  future  ships.  To  better  imderstand  ships'  responses  to 
wave  loading  trials  have  been  carried  out  by  DSTO  on  two  Royal  Australian  Navy 
vessels,  HMAS  SWAN  and  SYDNEY  so  as  to  monitor  and  record  ships'  motions  and 


1  Other  loading  types  such  as  grounding,  impact  and  blast  loading  are  not  relevant  to  this  discussion  and 
are  ignored. 
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structural  responses  in  a  seaway.  The  data  obtained  from  these  trials  is  used  for  a 
number  of  purposes  such  as  prediction  of  ship  responses  to  given  wave  loading 
conditions,  determination  of  the  effects  of  structural  modifications,  validation  of  finite 
element  models  and  to  provide  information  relating  to  the  fatigue  loading  of  RAN 
ships.  The  purpose  of  this  report  is  therefore  to  provide  an  outline  of  existing  analysis 
techniques  employed  and  the  manner  in  which  they  may  be  applied. 


2.  Wave  Analysis 

Ocean  waves  are  generated  primarily  through  the  action  of  the  wind  over  the  water, 
and  their  formation  is  influenced  by  such  factors  as  the  wind  strength  and  direction, 
the  time  for  which  the  wind  blows,  fetch  (the  distance  over  which  the  wind  blows)  and 
topography  of  the  ocean  floor.  The  waves  generated  form  an  irregular  pattern  of  wave 
heights,  lengths,  directions  and  frequencies  and  in  order  to  analyse  wave  records  and 
corresponding  responses  a  simplified  model  is  adopted.  Typically  spectral  analysis  is 
used  to  describe  the  wave  energy  as  a  function  of  wave  frequency  and  the 
corresponding  responses  at  these  frequencies  are  investigated  as  outlined  in  the 
following  sections. 


2.1  Encounter  Frequency 

As  a  ship  moves  through  a  seaway,  the  frequency  at  which  it  encounters  the  waves  is 
different  to  the  actual  wave  frequency  due  to  a  Doppler  shift  effect.  Figure  1  depicts  a 
ship  moving  with  velocity  at  an  angle  xj/  to  a  regular  seaway  with  circular  frequency 
to.  These  parameters  are  related  to  the  encounter  frequency  cOg  of  the  ship  by 

(Oj  =  0  -  (co^  /  g)Vs  cosxj/  .  (1) 
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The  encotinter  frequency  of  the  ship  is  then  dependent  not  only  on  the  wave  frequency, 
but  also  the  speed  of  the  ship  and  the  relative  angle  between  die  ship  and  the  seaway 
and  depending  on  these  parameters  a  number  of  alternative  situations  are  possible.  It 
can  be  seen  from  equation  (1)  that  in  bow  or  head  seas  ( 90°<\|f<270°)  cos(\i/)  is  always 
negative  so  the  encoimter  frequency  is  always  greater  than  the  wave  frequency.  If  \j/  = 
90°  or  \|r  =  270°  (beam  seas)  then  cos(\|/)  =  0  and  the  encounter  frequency  is  equal  to 
the  wave  frequency.  If,  on  the  other  hand,  the  ship  is  in  quartering  or  following  seas  (0 
°<\);<90°  or  270°<\|;<360°)  then  cos(\|r)  is  positive  and  so,  depending  on  the  wave 
frequency,  the  encounter  frequency  may  be  either  positive  or  negative. 

In  quartering  or  following  seas  the  encoimter  frequency  will  be  zero  if  co  =  0  or  co 
=g/ (Vg  cos\|/)  and  will  be  maximum  for  to  =g/(2Vg  cosxj; ),  in  which  case  =  g/(4Vs  cos 
Y  )  =  to/2.  In  deep  water  the  wave  velocity  (celerity)  c  =  g/to  so  that  high  frequency 
waves  advance  slowly  whereas  low  frequency  waves  advance  rapidly.  This  means 
that  if  the  encounter  frequency  is  zero  (to  0)  the  ship  will  be  travelling  at  the  same 
speed  as  the  waves,  if  encounter  frequency  is  positive  the  waves  will  be  overtaking  the 
ship,  and  if  encounter  frequency  is  negative  then  the  ship  is  overtaking  the  waves.  The 
relationship  between  actual  frequency  and  encounter  frequency  for  quartering  and 
following  seas  is  depicted  grapWcally  in  Figure  2.  It  is  seen  also  that  if  only  the 
absolute  value  of  encounter  frequency  is  known,  and  I  cOg  I  <  ojg  ,  as  might  be  the 
case  when  encounter  frequency  is  provided  by  on-board  instnunentation,  then  three 
wave  systems  are  possible. 
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CO 


Figure  2.  Encounter  Frequencies  in  Following  and  Quartering  Seas. 


Whatever  the  relative  direction  of  the  seaway,  the  possible  values  of  wave  frequency  co 
may  be  obtained  by  rearrangement  of  equation  (1),  which  gives 


CD  = 


g 

2Vg  cos\y 


1  ±V[1- 


4(D 
_ e 

g 


Vg  cosxy] 


rad  /  sec,  co  >  0 


(2) 


In  cases  where  more  than  one  solution  is  possible  (i.e.  more  than  one  real  solution)  it 
will  be  necessary  to  consider  other  information  such  as  on-board  observations  in  order 
to  determine  the  actual  value  of  co. 


2.2  The  Wave  Spectrum 

Waves  are  generated  through  the  action  of  the  wind  blowing  across  the  surface  of  the 
water.  Initially  ripples  are  formed  and  these  advance  with  increasing  fetch  until  waves 
are  created.  The  wind  will  continue  to  act  on  these  waves  generating  further  ripples 
and  waves  so  that  eventually  the  observed  waves  at  any  point  are  an  irregular 
superposition  of  waves  of  different  heights,  length  and  frequencies.  Under  ideal 
conditions  the  wave  components  all  travel  in  the  direction  of  the  wind,  but  this  is  never 
actually  the  case  and  it  is  foimd  that  the  irregular  sea  has  components  in  many 
different  directions,  a  phenomenon  referred  to  as  wave  spreading. 

A  seaway  is  an  essentially  random  process  and  it  is  assumed  that  the  irregular  sea  at  a 
particular  point  may  be  described  mathematically  by  the  linear  superposition  of  a  large 
number  of  regular  sinusoidal  waves  of  different  amplitudes  and  frequency,  that  is  as  a 
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Fourier  series.  Typically  this  information  is  presented  as  a  uni-directional  wave 
spectrum  which  provides  the  wave  energy  density  (energy  per  unit  horizontal  surface 
area)  as  a  fimction  of  wave  frequency.  TTie  complete  name  for  such  a  spectrum  is  the 
'Wave  Energy  Density  Spectrum,  however  it  is  usually  just  referred  to  simply  as  the  wave 
spectrum,  ^metimes  the  term  point  spectrum  is  used  to  denote  that  the  spectrum  is 
observed  at  a  fixed  point,  without  consideration  of  the  direction  of  component  waves. 
To  enable  components  from  different  directions  to  be  included  a  directional  wave 
spectrum  may  be  used  where  the  wave  energy  is  assumed  to  be  spread  over  a  range  of 
directional  components.  The  spreading  is  usually  assumed  to  be  S)anmetric  about  the 
primary  direction  of  wave  travel  (see  section  2.4)  and  so  does  not  accurately  describe  a 
bi-directional  sea  state  where  there  are  two  main  directional  components.  Also,  the 
superposition  of  regular  waves  does  not  exactly  describe  an  irregular  wave  system, 
particularly  in  higher  sea  states  where  non-linear  effects  such  as  wave  breaking  are 
more  pronounced,  however  it  is  relatively  simple  to  apply  and  for  many  applications  is 
quite  accurate. 

If  an  actual  ocean  wave  system  is  measured  as  a  discrete  time  record  {xr}  of  wave 
height  over  a  period  of  time  T  at  discrete  time  intervals  of  At,  the  actual  time  history 
may  be  approximated  as  the  sum  of  a  large  number  of  regular  sine  waves  at  discrete 
frequency  intervals  of  5(0  =  27i:  /  T.  A  Fourier  series  is  obtained  from  the  discrete  time 
series  using  the  Fast  Fourier  Transform  (FFT)  algorithm  where  the  Fourier  coefficients 
identify  the  magnitude  and  phase  lag  of  each  component  wave.  Although  the  ocean 
wave  system  is  considered  statistically  stationary  at  the  time  of  the  measurement,  the 
particular  time  history  being  examined  is  unique  and  will  never  be  repeated.  It  is  more 
useful  therefore  to  consider  the  energy  spectrum  obtained  from  the  Fourier 
representation  as  it  is  ideally  identical  for  adjacent  measurements  whereas  the  Fourier 
coefficients  of  the  component  waves  will  vary  from  sample  to  sample  due  to  random 
phasing.  The  energy  density  spectrum  is  then  obtained  essentially  by  forming  the 
complex  modulus  squared  of  the  Fourier  combinations,  although  additional  numerical 
manipulation  known  as  smoothing  is  required  to  provide  an  accurate  estimate  of  the 
true  spectrum.  The  steps  involved  are  explained  in  detail  and  in  summary  by  Newland 
(1984)  and  so  will  not  be  repeated  here. 

The  wave  spectrum  may  be  expressed  in  terms  of  either  circular  frequency  co  (rad/sec) 
or  cyclic  frequency  f  (hertz).  In  either  case  the  area  under  the  curve  must  remain 
unchanged  so  if  S^((0)  is  the  wave  spectrum  in  terms  of  circular  frequency  and 

if)  is  the  spectrum  in  terms  of  cyclic  frequency  then 

W^if  =  -^)  =2n.S^((0).  (3) 

A  typical  measured  wave  energy  density  spectrum  in  terms  of  cyclic  frequency  is 
shown  in  Figure  3.  This  spectrum  has  been  obtained  through  averaging  the  results  of 
20  overlapping  spectra,  each  obtained  using  2048  data  points. 
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Figure  3.  Typical  Measured  Wave  Energy  Density  Spectrum 

The  area  beneath  the  curve  of  an  infinitesimal  strip  of  width  6(0  (for  cyclic  frequency 
representation)  and  centred  at  is  proportional  to  the  energy  density  of  a  regular 
wave  of  frequency  (Bj^  .  This  proportionality  is  true  whichever  spectral  form  is  used. 
The  imits  of  the  spectral  density  ordinates  are  m2/rad/sec  and  since  the  energy  density 
of  a  regular  wave  of  frequency  co^  and  amplitude^  is  constant  of 

proportionality  must  therefore  have  units  of  pg  (kg/ m^  *m/ s^)-  The  area  beneath  the 
wave  spectrum  for  the  frequency  band  centred  at  co^  will  then  be  given  by^ 

pgS^((0k).6(0  =  ypgCo.k 

S^(o)k).5co  =  i  Co.k  m^/ rad /sec.  (4) 

The  right-hand-side  (RHS)  of  equation  (4)  is  in  fact  equal  to  the  variance  of  the  wave 
elevation  of  the  regular  sine  wave  of  frequency  coj.  and  amplitude  ,, .  That  is 


2  ^  is  used  to  denote  Wave  Elevation  while  denotes  Wave  Amplitude 

3  Care  should  be  taken  when  dealing  with  published  wave  spectral  data  as  in  some  cases  the  constant  of 


proportionality  may  be  taken  to  be  y  pg  in  which  case  S2^((0k;)  “  2S^((0|^)  — 


Si 

5co 
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E[Ck(t)']=S5((0,).5(0= 

The  total  variance  of  the  spectrum  is  obtained  by  integration  of  the  left-hand-side  of 
equation  (4).  Each  wave  component  may  be  assumed  to  be  independent  and  therefore 
provided  the  number  of  frequency  bands  is  sufficiently  large  so  that  8®  0,  the 

variance  of  the  sum  approaches  the  sum  of  the  variances  of  the  individual  component 
waves,  and 


E[C(t)"]=mo=|”s^(ro).d®.  (6) 

That  is,  the  area  beneath  the  spectrum  is  equal  to  the  variance  of  the  wave  system.  The 
notation  hiq  is  used  to  refer  to  the  variance  as  it  is  in  fact  the  zeroth  moment  of  the 
spectral  area  about  the  origin.  Since  the  variance  of  the  wave  elevation  is  proportional 
to  the  energy  density  of  the  wave,  mg  provides  a  good  measure  of  the  severity  of  the 
wave  system,  and  Vm^  is  in  fact  the  modal  or  most  probable  wave  amplitude  .  This 

is  not  to  be  confused  with  the  average  wave  amplitude  which  is  equal  to  1.25VmQ 
(see  further  in  section  2.5). 

If  the  spectrum  is  to  be  expressed  in  terms  of  the  cyclic  frequency  rather  than  the 
circular  frequency,  then 


and  again 

E[C(t)"]  W^CAd/  (7) 


2.3  Transformation  of  the  Wave  Spectrum 

The  wave  energy  spectrum  is  defined  for  a  fixed  point  in  space  and  so  it  is  frequently 
necessary  to  transform  the  spectrum  to  the  frame  of  reference  of  the  ship.  The  wave 
energy  spectrum  must  be  shifted  along  the  frequency  axis  so  that  the  frequency 
interval  of  8®  at  is  transformed  into  the  corresponding  encounter  frequency  interval 
8®g .  The  area  under  the  wave  energy  spectrum  must  be  the  same  in  each  case  so  that 
the  area  of  each  corresponding  incremental  strip  must  be  the  same.  That  is 

S5(®e,k).8®e  =  S^(®t).8® 


and  as  5®  and  8®^  become  infinitesimal 

S;(0)J  =  S^(®).d®/d®,. 
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So  from  equation  (1)  the  encounter  frequency  spectrum  is 


Sr  (©e)  =  Sr  (to). - - -  mVrad/sec.  (8) 

^  ^  g-2coVsCos\j/ 

Sometimes  it  is  convenient  to  consider  the  wave  system  in  terms  of  the  wave  slope 
spectrum,  for  example  when  considering  angular  motions  of  roll  pitch  and  yaw.  In 
deep  water  the  wave  slope  amplitude  is  given  by 

ao=©^Co/g  radians. 


And  in  the  same  way  as  the  wave  energy  spectrum  represents  the  energy  of  the  regular 
sine  waves  that  make  up  the  irregular  wave,  the  slope  spectrum  represents  the  slopes 
of  the  sine  waves  that  make  up  the  irregular  wave.  Then 


Sa(®k) 


2 


a 


1  0.k 

^  6o) 


S;(©k). 


(9) 


That  is,  the  wave  slope  spectrum  may  be  obtained  by  multiplying  the  wave  energy 
spectral  ordinates  by  co^/g^. 


2.4  The  Directional  Wave  Spectrum 

So  far  it  has  been  assumed  that  the  wave  energy  spectrum  has  been  derived  from  a 
record  of  the  surface  elevation  at  a  particular  point  without  regard  to  the  direction  of 
the  components  contributing  to  the  wave  system.  Invariably  the  wave  system  at  a 
particular  point  will  contain  contributions  from  several  directions,  and  the  amount  of 
spreading  at  any  particular  time  and  place  is  dependent  upon  both  the  topology  of  the 
ocean  floor  and  the  past  history  of  ^e  wave  system.  The  directional  wave  spectrum 
represents  both  the  frequency  and  the  directional  spread  of  energy  of  the  wave  system. 
It  is  defined  such  that  the  wave  energy  density  in  the  frequency  band  5(0  and  the 
directional  band  8(j)  is  equal  to  the  energy  density  of  the  regular  wave  of  frequency  cOj,. 
which  is  travelling  in  direction  (j); .  That  is,  similar  to  equation  (4), 

S5(©k.<l>)-S©.5<t)  =  ^  Co.k,i  mV  rad/sec.  (10) 

For  design  purposes  it  is  usual  to  assume  that  the  spreading  is  s)mimetric  about  the 
primary  wave  direction  and  is  distributed  over  a  spreading  range  of  ±<t)inax- 
achieved  through  the  use  of  a  cosine  spreading  factor  C((j))  such  that 
Sc((0^,(j))  =  C((|)).Sc((0t)  and  where  (Lloyd,  1989) 
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c((|))  =  D  cos^c— 

2<l>max 

=  0  elsewhere 

where  m  is  a  positive  integer  and 


■^max  <  <  ^max 


(11) 


D  = 


71 


4<t>a 


n 


1.3.5 . m _ 

2.4.6.  (m-1)  4(1)^ 
2.4.6.. ..m  7C 
1.3.5.  (m-1)  4^^ 


for  m  odd 


for  m  even. 


As  the  value  m  increases  the  energy  becomes  progressively  more  concentrated  around 
the  primary  wave  direction,  and  for  design  purposes  m  =  2  and  (|)^^  =  7c/2  are  usually 
chosen,  as  these  have  been  foimd  by  Cummins  and  Bales  to  be  appropriate  for  typical 
open  ocean  conditions  (Lloyd,  1989).  The  International  Ship  Structures  Congress  (ISSC) 
and  the  International  Towing  Tank  Conference  (l  ll'C)  both  recommended  values  of 
‘t>max  =  7t/2  for  maximiun  spreading  angle  but  with  m=4  and  m=2  respectively. 

The  total  wave  energy  distributed  over  the  range  ±(l)niax  point  spectrum  is  then 
obtained  by  integration  over  that  range. 

S,(a))  =f7  S,(co,(|)).5(l)  (12) 


In  practice  equation  (12)  is  evaluated  at  discrete  intervals  of  say  15°  and  then  the  point 
spectrum  is  obtained  by  numerical  integration  over  the  spreading  range. 


2.5  The  Statistics  of  Irregular  Waves 

Over  a  period  of  20-30  minutes  a  wave  system  at  a  particular  point  may  be  considered 
statistically  stationary  and  so  certain  statistical  characteristics  may  be  determined. 
Wave  parameters  commonly  of  interest  are: 

•  wave  amplitude  (metres)  and  the  mean  value  of  many  measurements 

•  wave  height  H  (metres)  and  the  mean  value  of  many  measmements  H 

•  wave  period  Tp  between  adjacent  peaks  (seconds)  and  the  mean  value  of  many 
measurements  Tp 
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•  zero  crossing  period  Tz  (seconds)  and  the  mean  value  of  many  measurements 
Tz 

Wave  height  and  amplitude  are  usually  defined  as  being  for  each  event  ^tween  zero 
up-crossings  and  it  is  also  common  to  refer  to  a  significant  wave  height  H1/3  which  is 
the  mean  value  of  the  l/3rd  highest  of  many  wave  height  measurements.  A  typical 
fixed-point  wave  record  and  the  various  wave  parameters  are  shown  in  Figure  4. 


Maximum 

Trough 


Figure  4.  Typical  Fixed-Point  Wave  Record  and  Definition  of  Parameters 

In  equation  (6)  the  variance  was  obtained  by  taking  moments  of  the  wave  spectrum 
about  the  origin  to  =  0.  In  a  similar  manner,  higher  moments  of  the  wave  spectrum 
may  be  determined  and  used  to  compute  a  number  of  wave  parameters  as  follows 
(Bishop  and  Price,  1979). 

mg  =  variance  of  the  wave  elevation 

=  zeroth  moment  of  the  wave  spectrum 

=  S^(o)).d0  (m^) 

=  the  area  tmder  the  wave  energy  density  spectrum 
=  mean  square  wave  elevation  (for  zero  mean  process) 

m2  =  second  moment  of  the  wave  spectrum 

=  S^(a)).dco  =  o)^S^(co).da)  (mVsec^) 

=  the  area  xmder  the  wave  vertical  velocity  spectrum 
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m4  =  fourth  moment  of  the  wave  spectrum 

=  J“  Sj(co).da)  =  (D‘^S5(o)).dco  (mZ/sec^) 

=  the  area  under  the  wave  vertical  acceleration  spectrum 

lUj^  =  n'th  moment  of  the  wave  spectrum 

=  co”S^(co).dto  (m2/sec") 

0)  =  average  frequency 

=  mj  /mo  (rad/sec) 

T  =  characteristic  period  (of  the  component  wave  of  average  frequency  CO ) 

=  27C  mo/mi  (sec) 

Tp  =  average  wave  period  between  peaks  or  hollows 

=  27rV(m2/m4)  (sec) 


Tz  =  average  zero-upcrossing  period 

=  27cV(mo/m2) 

Cm  =  most  probable  (modal)  wave  amplitude 

=  ^1% 

Ca  =  average  wave  amplitude 

=  1.25Vmo 


H  = 

mean  peak  to  trough  wave  height 

2.5lVmo 

(m) 

Hi/3  = 

significant  wave  height  (average  of  l/3rd  highest  waves) 
4.0Vmo 

(m) 

Hmo  = 

average  height  of  1/lOth  highest  waves 

S.OSVmo 

(m) 

1  = 

average  number  of  zero  upcrossings  per  emit  time 

^VfiHo/mj) 

271 

The  formulae  for  average  periods,  heights  and  amplitudes  are  strictly  only  valid  if  the 
wave  elevation,  measured  at  equal  time  intervals,  is  a  zero-mean  process  and  has  a 
Normal  or  Gaussian  distribution,  and  the  wave  spectrum  is  narrow  banded.  This 


(sec) 

(m) 

(m) 
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occurs  if  the  phase  angles  of  the  component  regular  waves  have  uniform  probability 
between  0  and  2jr,  and  fortunately  in  the  short  term  (possibly  up  to  several  hours  in 
steady  conditions),  this  is  approximately  true  for  ocean  waves. 

The  bandwidth  parameter  is  defined  as 


and  is  usually  small  and  any  errors  associated  with  the  above  formulas  for  wave 
height,  period,  etc.  are  also  quite  small.  Ochi  (1981)  showed  for  a  large  number  of  wave 
cycles  that  when  determining  extreme  values  with  a  low  probability  of  exceedance 
(section  5.3),  the  error  associated  with  this  assumption  is  quite  small  even  for  values  of 
e  up  to  0.9. 

Empirical  relationships  have  been  developed  which  relate  significant  wave  height  and 
zero-upcrossing  period  with  visual  estimates  of  wave  heights  and  periods.  It  is 
generally  considered  that  visually  estimated  wave  heights  can  be  taken  as  the 
significant  wave  height,  for  a  long  time  it  was  generally  accepted  that  observed  period 
corresponds  to  the  zero  crossing  period.  Visual  estimates  of  wave  period  are  generally 
considered  to  be  inaccurate,  particularly  if  the  observed  sea  is  bi-directional  and  made 
up  of  a  combination  of  swell  and  wind  waves.  The  following  relationships  were 
derived  by  Nordenstrom  from  a  large  number  of  visual  observations  and  are 
sometimes  used  to  improve  the  accuracy  of  visual  observations  (Lewis,  1989).  It  should 
of  course  always  be  remembered  that  there  is  a  considerable  scatter  in  the  data  used  in 
their  derivation,  and  that  visual  wave  estimation  is  a  subjective  assessment. 

Hi/3  =  1.68 

%  =  0.82  T^-^® 

(Hy  andTy  are  the  visual  estimates  of  wave  height  and  period  respectively.)  More 
recently  it  has  been  found  (Hogben  et  al,  1986)  that  visually  observed  wave  period 
corresponds  more  closely  to  the  modal  wave  period,  Tp  ( the  period  corresponding  to 
the  peak  of  the  wave  energy  spectrum),  and  that  it  is  related  to  zero-crossing  period  by 
Tp  =  1.4  Tz. 

Joint  probability  distributions,  or  scatter  diagrams,  of  wave  heights  and  wave  periods 
for  104  areas  covering  the  globe  are  provided  by  Hogben  (Hogben  et  al,  1986).  The 
results  of  55  million  observations  are  sxunmarised  and  are  enhanced  to  accoimt  for 
wind  speed  and  direction.  The  scatter  diagrams  provide  the  joint  probabilities  of  wave 
heights  in  increments  of  one  metre  and  wave  periods  in  increments  of  one  second  and 
are  broken  down  into  geographical  area,  season  and  45°  directional  groupings. 
Seasonal  groupings  depend  on  the  local  seasonal  patterns  and  directional  groupings 
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define  the  direction  of  the  wind  generated  wave  component.  Wave  height  and  period 
given  in  the  tables  correspond  to  the  significant  wave  height  and  zero-crossing  period 
respectively.  When  processing  the  raw  visual  data  only  those  wave  height 
observations  which  included  both  wind  waves  and  swell  were  included  and  these 
were  combined  to  give  the  resultant  significant  wave  height  Hi/3=(H^  + 

Visual  observations  of  wave  period  were  not  used  as  these  are  considered  to  be 
unreliable.  Instead  wave  period  statistics  based  on  the  modelling  of  joint  probability  of 
wave  period  with  wave  height  derived  from  the  analysis  of  instrument  measured  data 
are  used. 


3.  Idealised  Wave  Spectral  Families 

Various  attempts  have  been  made  to  establish  idealised  wave  spectral  families  which 
would  enable  designers  and  hydrographers  to  define  the  wave  spectrum  through  a 
limited  number  of  parameters.  Often  the  use  of  these  spectra  is  appropriate  only  to  a 
particular  stage  of  wave  generation  so  care  must  be  taken  to  use  them  accordingly. 
Most  of  these  spectra  are  uni-modal  (singular  peak)  and  take  the  form 

=  ^exp[-B/0P]  (15) 

^  0)“ 

where  a  and  P  are  empirically  derived  coefficients  while  A  and  B  are  constants.  This 
may  also  be  written  in  a  non-dimensional  form,  relative  to  the  modal  value  of  the 
spectrum S^(co^),  the  maximum  value  given  by  dS^  (co)  /  dco  =  0  and  at  which  point 

co^j  =  (Bp  /  a)*^^.  Then  (Bishop  and  Price,1979) 

=  (co/(o,)-“expfc[l-(G)/a)„)-P]| 

IP  J 


3.1  Pierson-Mokowitz  Spectrum 

The  Pierson-Mokowitz  spectrum  is  an  open  water  spectrum  of  the  general  type  of 
equation  (15)  but  which  only  depends  upon  a  single  parameter  V^ ,  the  wind  speed.  It 
is  intended  to  describe  the  point  spectrum  for  fully  developed  seas  and  was  developed 
primarily  for  oceanographic  use.  It  is  perhaps  of  limited  use  to  the  naval  architect 
because  for  fully  developed  seas  to  occur  there  must  be  large  fetch,  constant  wind 
speed  and  direction  over  several  hours  and  no  contamination  from  swells  in  other 
wave  generation  areas.  The  spectrum  is  given  by 

S^(0))  =  ^exp[-B(g/V,)'‘]  (16) 

^  0) 
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where  S^(co)  =  the  spectral  ordinate  in  cm^.sec 

CO  =  frequency  in  radians/ sec 

A  =  8.10  *  10-3 

B  =  0.74 

=  wind  speed  in  cm/ sec  measured  at  19.5  m  above  the  sea  sirrface 

An  estimate  of  the  average  wave  period  (UCL,  1988)  using  this  spectrum  is  given  by 

T  =  0.8l(27tV„/g) 

and  of  significant  wave  height  by 

Hjy3  =  2.12  X 10’^ 

where  is  in  metres/sec. 

3.2  The  Bretschneider,  ISSC  and  ITTC  Spectra 

The  Bretschneider  spectrum  is  also  an  open  water  spectrum  of  the  kind  given  in 
equation  (15)  with  a  =  5  and  p  =  4,  and  is  designed  to  represent  both  rising  and  falling 
seas  as  well  as  fully  developed  seas.  While  individual  point  spectra  may  not  match  the 
Bretschneider  spectrum  particularly  well,  it  was  found  by  Ochi  and  Bales  that  the 
differences  tended  to  disappear  through  statistical  averaging  (Lewis,  1989).  It  is  at  best 
to  be  considered  approximate  but  it  does  permit  both  period  and  wave  height  to  be 
specified  separately.  The  spectrum  is: 

Sk(co)  =  exp[-B/co‘‘]  (17) 

and  the  two  constants  A  and  B  are  given  by 

B  =  1.25 
A  =  4mo  B 

This  spectrum  was  adopted  by  the  ISSC  (1967)  to  be  used  in  conjunction  with  observed 
wave  height  and  period,  which  are  taken  to  be  the  sigmficant  wave  height  and  the 
characteristic  period  as  given  in  section  2.5,  and  where 

A=173h)^/T^  (18) 

B  =  691/  f'* 
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The  characteristic  period  T  is  usually  taken  as  the  zero  crossing  period  Tr.  When 
visually  observed  wave  height  and  period  are  used  for  equation  (18)  it  is 
recommended  that  the  results  given  in  equation  (14)  be  used  to  determine  significant 
wave  height  and  zero  crossing  period  before  calculating  the  spectnun  according  to 
equations  (17)  and  (18). 

The  Bretschneider  spectrum  was  also  recommended  by  the  15th  International  Towing 
Tank  Conference  (ITTC,  1978)  for  average  conditions  (not  fully  developed  seas)  when 
the  significant  wave  height  is  the  only  information  available,  then  in  equation  (17) 


A  =  8.10  *  10’^ 
B  =  3.11/Hj^ 


(19) 


3.3  Joint  North  Sea  Wave  Project  Spectrum 

Joint  North  Sea  Wave  Project  (JONSWAP)  was  set  up  to  study  the  growth  of  waves 
under  conditions  of  limited  fetch  and  under  the  effects  of  shoaling  water.  A  family  of 
spectra  was  computed  from  the  measurements  made  by  the  project  and  are  again  of 
the  form  given  by  equation  (15)  but  modified  by  a  frequency  dependent  (and 
somewhat  complicated)  peak  enhancement  factor.  The  effect  of  the  peak  enhancement 
factor  is  to  produce  a  spectrum  which  is  narrower  and  has  higher  peaks  than  the 
standard  Pierson-Mokowitz  form,  and  which  is  more  suitable  for  storm  conditions 
than  fully  developed  open  ocean  waves.  The  spectrum  is  dependent  upon  wind  speed 
and  fetch,  and  although  it  can  at  times  provide  a  better  fit  to  data  than  other  spectra  it 
does  not  lend  itself  to  practical  use  for  design  or  most  response  analysis  purposes.  A 
modified  form  of  the  JONSWAP  spectrum  has  been  proposed  by  Det  Norske  Veritas 
(UCL,  1988)  which  enables  the  spectrum  to  be  calculated  using  significant  wave  height 
and  zero  crossing  period. 


8^(0)) 


h?/3t,  rcoT.y^ 

871^  I  27t  j 


7t 


(20) 


3.4  The  Ochi  6-Parameter  Spectrum 

A  family  of  spectra  were  developed  by  Ochi  and  Hubble  which  generally  provides  a 
much  better  fit  to  measmed  wave  spectra  than  do  the  previously  mentioned  spectra, 
particularly  in  cases  where  the  spectrum  is  bi-modal  (Lewis,  1989).  This  spectral  family 
has  the  form 
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S^Cco) 


( 


4^+1  \X 


nx) 


/co)^] 


0) 


4X+1 


(21) 


where  H,  is  the  significant  wave  height,  r(^)  is  the  gamma  function  and  X  is  a  shape 
parameter.  The  spectrum  given  by  equation  (21)  is  um-modal  with  the  peak  becoming 
sharper  with  increasing  X.  By  adding  two  of  these  forms  together  a  six  parameter 
spectrum  is  obtained  which  was  able  to  better  match  commonly  observed  bi-modal 
spectra.  Then 


1  '  A  TO.]' 

Si«»)  =  J  S  - 

^  j 


H 


V3,j 


ra) 


4Xi  +1 
CO  J 


4X ;  +  1  , 

■exp[ - - (co„j/®)'^] 


(22) 


where  j=l,2  stands  for  the  lower  and  higher  frequency  components,  respectively.  The 
six  parameters  H  1/3,1,  Hj/3.2,  CO^i,  0)^,  ^i  and  X^  may  be  determined  numerically  by 

minimising  the  difference  between  the  measured  spectrum  and  the  calculated 
spectrum.  Although  bi-modal  spectra  are  able  to  be  represented,  no  mformation  is 
given  as  to  the  relative  directiorvs  of  the  waves.  For  design  purposes  it  is  considered 
more  appropriate  to  use  the  ISSC  spectrum,  as  although  it  may  not  be  as  accurate  for 
specific  spectra,  it  adequately  describes  the  statistical  average  and  is  simpler  to  apply. 


4.  Response  Amplitude  Operators 

The  relationship  between  the  wave  spectrum  and  the  ship  response  (typically  strain  or 
motions)  spectrum,  835  (coj,  is  often  assumed  linear  and  is  expressed  in  terms  of 
Response  Amplitude  Operators.  The  RAO  is  a  frequency  dependent  function  which 
provides  the  ratio  of  the  response  spectrum  to  the  input  spectrum  over  the  entire 
frequency  range.  It  is  the  modulus  squared  of  the  frequency  response  function  or 
transfer  function,  and  does  not  therefore  contain  any  phase  information  between  the 
response  and  the  input.  As  the  ship's  response  is  always  in  terms  of  encoxmter 
frequency  the  relationship  is  written  as 

Sgi(0)J  =  IHgi^  (co,)P  .S^f©,)  =  RAO.S^ffflJ  (23) 

However  as  has  been  shown  in  Section  2.1  there  is  not  always  a  one-to-one  relationship 
between  the  encounter  and  actual  frequency.  For  the  range  of  headings  between  90°< 
\ii<270°  there  is  a  one-to-one  relationship  in  the  transformation  from  ©  to  tcfe,  but  for 
quartering  or  following  seas  this  is  not  necessarily  the  case.  It  can  be  seen  from 
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equation  (8)  that  as  co  approaches  g  /  2V5Cos\|/  the  denominator  approaches  zero  and 
hence  a  singularity  occurs.  The  evaluation  of  the  spectra  and  RAO  must  be  broken 
down  into  three  stages  as  depicted  in  Figure  2  and  a  full  description  of  the  process  is 
provided  by  Lewis. 

RAOs  may  be  calculated  either  theoretically  using  strip  theory  programs,  or  may  be 
derived  experimentally  by  monitoring  of  ships  or  ship  models.  The  concept  of  wave 
spreading  was  raised  in  section  2.4  where  the  use  of  a  spreading  factor  was  applied  to 
the  wave  spectrum  to  allow  for  the  different  directions  of  wave  components.  When 
evaluating  the  response  spectrum  the  directional  wave  spectrum  may  be  used  to  allow 
for  wave  spreading,  so  then 

SgtfCOe)  =  ((OjP  .S^((D,(l)).8<t)  .  (24) 

In  practice  equation  (24)  is  evaluated  at  discrete  intervals,  of  say  15°  and  the  total 
response  spectrum  is  obtained  by  numerical  integration  over  the  spreading  range. 

As  mentioned  in  section  2.2  the  use  of  a  linear  wave  spectrum  does  not  exactly  describe 
an  irregular  wave  system,  particularly  in  higher  sea  states.  Similarly  the  response  of 
a  ship  to  the  seaway  is  not  exactly  linear,  especially  to  the  higher  sea  states,  and  so  the 
response  spectrum  provided  in  equation  (24)  will  not  be  exactly  correct.  When  strip 
theory  programs  are  used  during  the  design  process,  the  calculation  of  the  RAOs 
involves  even  more  approximations  as  further  simplifying  assumptions  are  used  in  the 
calculations  of  the  wave  spectrum  and  the  RAOs.  Furthermore  strip  theory 
calculations  do  not  usually  differentiate  between  hogging  or  sagging  responses  which 
are  in  practice  quite  different.  The  end  result  is  that  tiie  calculated  response  spectrum 
typically  over  estimates  the  response  of  the  ship  for  higher  wave  loadings,  particularly 
the  hogging  response  (Clarke  1986).  One  approach  to  try  to  improve  these  predictions 
is  to  carry  out  either  model  or  full  scale  tests  on  ships  to  monitor  their  response  to  the 
seaway  and  thereby  determine  experimental  RAOs.  These  RAOs  may  then  be  used 
instead  of  those  from  strip  theory  calculations  to  determine  the  ship's  response  to  a 
particular  seaway,  or  alternately  to  modify  the  calculated  RAOs  so  that  the  strip  theory 
program  may  more  accurately  determine  the  ship's  response.  This  modification 
would  typically  take  the  form  of  a  scaling  factor,  and  it  would  be  possible  to  use  this 
scaling  factor  in  the  assessment  of  responses  of  geometrically  similar  ships. 


5.  Probability  Formulae 

So  far,  the  discussion  has  revolved  aroimd  only  short  term  assessment  of  sea  state 
conditions  and  the  corresponding  ship's  responses.  As  the  sea  state  as  well  as  the 
ship's  speed,  position,  orientation  and  other  loading  conditions,  are  continuously 
changing,  it  is  impossible  to  know  in  advance  with  any  certainty  all  conditions  that  the 
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ship  will  encoimter,  and  even  if  it  were,  it  would  be  totally  impractical  to  assess  ^e 
ship's  response  to  all  of  these  conditions.  To  be  able  to  design  safe  shructures  e 
designer  is  therefore  required  to  make  predictions  regarding  the  lifetime  loading 
conditions  such  as  the  fatigue  cycling  and  extreme  loads  on  the  ship.  As  wave 
loading  is  of  a  random  nature  these  predictions  are  based  upon  wave  and  respoi^e 
statistics  and  so  an  imderstanding  of  the  relevant  probability  theory  and  its 
application  is  required. 


5.1  Short  Term  Probabilities 

As  already  mentioned  in  section  2.5  the  short  term  distribution  of  wave  elevation  can 
be  approximately  described  by  the  Gaussian  distribution  provided  the  spectrum 
bandwidth  is  narrow.  Likewise  the  linear  response  of  the  ship,  such  as  strain  or  ship 
motions,  can  be  assumed  to  follow  a  Gaussian  distribution.  Figure  5  shows  a  typical 
histogram  of  strain  response  on  the  main  deck  of  a  frigate  measured  at  equal  time 
intervals  over  a  period  of  several  minutes.  The  corresponding  Gaussian  distribution 
curve  is  also  plotted  and  it  is  seen  that  while  the  curve  provides  a  reasonable 
approximation  to  the  data,  it  is  not  a  true  Gaussian  process.  As  is  quite  often  the  case 
the  level  of  sagging  responses  is  higher  than  for  the  hogging  responses  due  to  non¬ 
linear  effects  such  as  flared  bow  and  stem  sections.  The  measured  response  may  also 
be  affected  by  the  location  of  the  strain  gauge  as  shear  lag  effects  may  reduce  the 
apparent  compression  loading. 


0 


Strain  Response 

Figure  5.  Strain  Response  Histogram  and  Corresponding  Normal  Distribution 
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The  Gaussian  probability  density  function  of  a  continuous  random  process  is  given  by 


f(x)  = 


1  g-(x-p)^/2o^ 

o-s/^ 


(25) 


where  p  is  the  mean  and  a  is  the  standard  deviation  of  the  process.  In  relation  to  wave 
elevation  and  corresponding  response  analysis  the  process  will  have  a  zero  mean 
valued  so  then  the  probability  that  an  individual  measurement  ^  =  x/a  will  fall 
between  the  values  4i  and  ^2  is  given  by 


Intergration  of  the  RHS  of  the  equation  is  difficult  and  so  values  are  obtained  from 
tables. 

Usually  it  is  the  peak  value  of  the  process  rather  than  the  variation  with  time  that  is  of 
interest,  and  whereas  the  short  term  distribution  is  approximated  by  the  Gaussian 
distribution,  it  is  found  that  the  short  term  distribution  of  the  peak-to-trough  (double) 
amplitudes  is  approximated  by  a  Rayleigh  distribution. 

fft.)  =  (26) 


Note  that  in  relation  to  a  ship's  responses  to  ocean  waves,  there  will  actually  be  a 
different  distribution  for  the  positive  and  the  negative  amplitudes  if  they  are 
considered  individually.  In  other  words  the  hogging  and  sagging  responses  are 
different,  and  this  is  well  documented  (e.g.  Clarke,  1986).  However  for  the  present  it 
is  assumed  that  the  Rayleigh  distribution  is  of  the  peak-to-trough  amplitude,  or  double 
amplitude,  so  that  the  mean  value  of  the  process  is  zero.  If  the  single  amplitude 
distribution  is  to  be  obtained  it  must  be  remembered  that  the  mean  value  is  no  longer 
zero.  The  probability  that  the  double  amplitude  will  fall  between  and  ^2  is  then 


^  It  is  assumed  that  in  the  analysis  of  response  data,  such  as  recorded  strain  gauge  data,  the  analysis 
begins  by  adjusting  the  data  values  so  that  a  zero  mean  process  is  achieved. 
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while  the  probability  that  the  double  amplitude  will  exceed  a  certain  value  gives  the 
short  term  probability  of  exceedance  in  one  cycle 


Psi  =  F(^a>y 


(27) 


5.2  Long  Term  Probabilities 

What  is  usually  of  interest  is  the  ship's  long  term  response  to  loading  over  a  wide 
range  of  operating  conditions  and  environments,  so  as  to  be  able  to  predict  the  fatigue 
cycling  and  expected  maximum  lifetime  loading  on  the  vessel.  The  long  term 
cumulative  probability  distributions  of  loading  is  therefore  required.  To  achieve  this  it 
is  first  necessary  to  obtain  and  analyse  a  large  database  of  information  relating  to 
various  sea  states  likely  to  be  encountered,  as  well  as  information  relating  to  the  ship  s 
responses  to  these  sea  states.  For  this  purpose  monitoring  trials  are  often  carried  out 
on  ships,  (e.g.  see  appendix  2  of  Lewis  and  Zubaly,  1981)  and  then  the  data  must  be 
analysed  so  as  to  determine  the  long  term  probability  distribution.  Typically  the 
following  steps  will  be  taken: 

(i)  Full  scale  ship's  trials  are  conducted  over  a  broad  range  of  sea  state  and  ship's 
operating  conditions  to  obtain  a  comprehensive  set  of  short  term  response  data. 

(ii)  Each  short  term  record  is  summarised  using  histograms  of  response  level 
versus  frequency,  and  the  spectral  properties  and  Rayleigh  distribution  for  each 
may  be  calculated. 

(iii)  The  cumulative  probability  distributions  are  then  determined  based  upon 
either  the  modal  (most  probable)  or  extreme  wave  height  values  from  these 
records.  Commonly  a  theoretical  distribution,  such  as  the  Weibull  distribution, 
is  fitted  to  the  calculated  distribution  for  the  purpose  of  extrapolation  to  longer 
term  responses. 


5.2.1  Histogram  Analysis 

Histograms  enable  the  large  amount  of  data  collected  during  the  short  term  trials  to  be 
summarised  by  recording  the  measured  responses  in  certain  pre-determined  ranges 
against  their  frequency  of  occurrence.  By  measuring  response  cycles  over  a  long 
period  of  time  and  a  range  of  operating  conditions,  a  history  of  responses  and  their 
relative  frequency  of  occurrence  is  obtained.  The  histograms  provide  the  probability 
density  of  occurrence  of  particular  response  levels  while  the  cumulative  distribution 


5  It  is  to  be  remembered  that  the  probability  density  function  may  be  obtained  by  differentiation 
of  the  cumulative  probability  distribution. 
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can  be  obtained  by  numerical  integration  of  the  histogram  data.  If  all  of  the  histogram 
data  from  all  of  the  trials  is  combined  the  cumulative  distribution  gives  a  plot  of 
response  level  versus  probability  of  occurrence  over  the  life  of  the  trials.  When  plotted 
on  a  semi-logarithmic  plot  the  data  will  often  approximately  follow  a  straight  line, 
except  possibly  at  the  lowest  response  levels.  The  distribution  obtained  through  this 
approach  is  limited  to  the  time  over  which  the  data  is  collected  and  should  not  be  used 
directly  for  extrapolation  to  longer  periods  as  the  results  are  limited  to  the  actual 
weather  and  sea  conditions  encountered,  the  ship's  structural  design  and  characteristic 
response  to  the  seaway.  The  results  obtained  are  therefore  biased  towards  the  actual 
conditions  of  the  trials  and  do  not  necessarily  give  correct  weighting  to  all  sea  and 
operating  conditions  which  may  be  encoxmtered  during  the  life  of  the  ship.  Also, 
comparative  trials  have  shown  that  responses  vary  between  similar  ships  of  the  same 
class  and  even  between  port  and  starboard  locations  on  the  same  ship  (e.g.  Hoffman 
and  Lewis,  1969). 

The  histograms  are  therefore  useful  in  the  analysis  of  existing  ships  with  long  histories 
of  voyage  monitoring  but  should  not  be  applied  directly  to  other  ships  or  in  ship 
design  unless  variations  in  probable  sea  and  operating  conditions  can  be  properly 
accoimted  for.  It  may  also  be  necessary  to  make  certain  assumptions  as  to  the  form  of 
the  long  term  distribution  so  as  to  be  able  to  make  consistent  predictions  of  the  long 
term  responses  and  two  common  methods  are  used;  extrapolation  based  upon  either 
modal  or  extreme  wave  height  values  (Liu  et  al,  1981). 


5.2.2  Extrapolation  Based  upon  Modal  Values 

The  extrapolation  from  measured  responses  to  longer  periods  may  be  based  upon  the 
assumption  that  the  short  term  distribution  of  peak  values  within  each  record  follow 
the  Rayleigh  distribution  and  that  the  modal  or  most  probable  values  of  the  records 
themselves  follow  some  form  of  probability  distribution.  For  a  zero-mean  process  the 
modal  amplitude  is  equal  to  the  standard  deviation  a  and  so  completely  defines  the 
Rayleigh  distribution.  More  typically  the  modal  or  most  probable  wave  height  value  is 
used  when  describing  the  sea  conditions  and  so  the  long  term  probability  density 
function  of  the  modal  wave  height  values  may  be  given  by  g(-\^H  )  where 

■JEh  =  72^.  The  long  term  probability  density  function  is  defined  to  accovmt  for 
the  variations  due  to  different  sea  states,  the  ship 's  speed,  heading,  loading  condition 
and  so  on^. 

The  joint  long  term  probability  p(x,Ve  h  )  is  the  joint  probability  of  a  particular  value  x 
(typically  bending  moment  or  stress  level)  with  VEh  ,  and  is  by  definition 


^  The  overscore  ~  is  used  to  indicate  that  g  (VEh  )  is  the  probability  distribution  obtained  from  the 
recorded  data. 
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p(x,7^)  = 

where  the  probability  density  function  f(x|.^^^ )  is  the  conditional  probability  of  x  for 
given  .Je^  and  is  assumed  to  follow  the  Rayleigh  distribution. 

Different  theoretical  forms  have  been  suggested  for  gCVEn  )  such  as  the  Gaussian, 
Weibull  and  Gumbel  distributions.  There  appears  little  theoretical  justification  for  ^y 
of  these  forms,  other  than  that  they  work  reasonably  well,  although  the  Gaixssian 
distribution  tends  to  over  estimate  the  larger,  less  frequent  values,  while  the  Weibull 
distribution  gives  a  greater  weighting  to  the  lower  values.  The  Weibull  distribution 
appears  to  be  gaining  broad  approval  as  it  is  found  that  much  of  the  fatigue  damage  on 
ships  is  caused  by  these  more  moderate  loads  which  occur  with  greatest  regularity. 
However  it  should  be  noted  that  it  is  not  essential  that  gC-^E^j  )  be  fitted  to  any 
particular  theoretical  distribution. 

The  long  term  probability  that  x  will  exceed  the  level  is  then  given  by  substituting 
the  Rayleigh  function  into  equation  (28)  and  integrating  with  respect  to  x,  which  gives 


If  gfVEH  )  is  assumed  to  follow  one  of  the  above  standard  distributions,  certain 
parameters  which  define  the  distribution  must  be  evaluated  and  this  is  done  by  fitting 
the  recorded  data  to  the  selected  distribution  using  a  suitable  error  minimisation 
technique.  Numerical  evaluation  of  equation  (29)  is  then  usually  required,  for  which 
an  upper  limit  of  five  to  six  times  the  standard  deviation  has  been  shown  to  provide 
the  minimum  value  to  ensure  sufficient  accuracy  to  the  final  result  (Hoffman  and 
Lewis,  1969). 

Evaluation  of  the  integral  (29)  provides  the  long  term  probability  that  on  average  the 
level  Xq  wiU  be  exceeded  just  once  in  the  lifetime  of  the  ship  as  illustrated  by  point  A  in 
Figure  6.  But  if  many  similar  ships  are  considered  it  would  be  foimd  that  in  some 
cases  the  level  Xp  will  not  be  exceed  while  in  other  cases  Xq  will  be  exceed  more  than 
once.  In  fact  there  is  a  67%  (Lewis  and  Zubaly,  1981)  probability  that  this  will  be  the 
case,  which  is  unacceptably  high  for  design  purposes.  A  probability  that  the  given 
level  Xp  is  not  exceeded  except  with  an  acceptably  low  level  of  risk  (such  as  0.01)  by 
any  one  of  a  large  number  of  ships  of  the  same  class  is  therefore  required.  This  design 
lifetime  probability,  Pln  is  given  by  the  product  of  the  cumulative  probability  at  the 
design  level  x^  times  the  number  of  cycles  and  illustrated  by  point  B  in  figure  6.  The 
level  Xp  will  then  not  be  exceeded  over  the  lifetime  of  the  vessel  with  a  probability  of 
Pjjj  =  N*Pli(xd).  For  example  if  a  value  such  as  N=3.0x10'7  is  chosen  with  a  lifetime 
design  probability  of  0.01,  then  the  design  response  value  x^  is  obtained  reading  off 
the  curve  at  Pli(x)  =  Pln  /N  i.e.  at  Pu  =  3.0x10"^ 
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No  of  Cycles  N 


Figure  6.  Long  Term  Probability  of  Exceedance  Based  Upon  Modal  Values 


5.2.3  Extrapolation  Based  upon  Extreme  Values 

Typical  of  common  engineering  practice,  extrapolation  to  determine  lifetime  loading 
may  also  be  based  upon  extreme,  rather  than  modal,  values.  This  approach  requires 
potentially  less  data  collection  effort  as  it  is  only  the  extreme  value  from  each  data 
sample  which  is  of  interest.  However  the  method  may  also  be  less  reliable  as  not  all  the 
sample  information  is  used,  although  that  this  is  so  is  not  clearly  the  case. 

If  the  number  of  data  samples  from  short  term  trials  is  n  then  a  probability  distribution 
of  the  extreme  values  may  be  determined.  For  each  sample  there  is  a  probability 
F,(Xi  <  Xq)  that  the  extreme  response  Xj  is  less  than  a  particular  value  Xq  in  any  one 
cycle.  So  for  the  n  samples  the  probability  that  each  of  these  extreme  response  in  any 
one  cycle  is  less  than  Xq  is  the  same  as  the  probability  that  the  maximum  value 
Ti=  max(Xi,X2,...X„)  is  less  than  Xq  and  is  given  by 

G(ti<Xo)  =  Fi(X,  <Xo)nF2(X2  <Xo)n../^F,(X„  <Xo)  (30) 

That  is,  the  cumulative  probability  distribution  G(x)  is  given  by 
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G(x<Xo)=  F“(Xo)  (31) 

where  ^  (xq)  is  a  shorthand  way  to  denote  the  combined  probability  distribution 
functions  of  the  individual  variables,  each  of  which  are  evaluated  at  the  chosen  peak 
value  Xq.  Similarly  the  probability  that  each  extreme  response  is  greater  than  Xq  will 
be 


G(x>Xo)  =  [l-^fxo)]"  .  (32) 

The  result  given  in  (32)  is  termed  the  long  term  probability  of  exceedance  in  one  cycle  and  is 
denoted  Pli-  What  is  needed,  however,  is  the  probability  that  x  is  greater  than  Xq  at 
least  once  in  a  lifetime  of  N  cycles.  So,  assuming  n  is  large  enough  to  enable  the  long 
term  distribution  G(x)  to  be  given  by  F”(Xo),  by  defining  the  probabiUty  per  cycle 
G(x>Xo)  =  Pli  as  "success"  and  G(x<Xo)  =  l-Pu  as  "failure"  then  by  application  of 
the  Binomial  distribution  it  can  be  shown  that  for  large  N  that  the  long  term  probability 
of  exceedance  in  N  cycles  is 


Pm(x>Xo)  =  1  -  .  (33) 

It  can  be  seen  from  equation  (33)  that  if  Pn  =  1/N  then  Pj^,  =  1  -  e  ^  =  0.632,  so  that 
in  a  similar  fashion  to  section  5.2.2.  there  is  a  63%  chance  that  the  value  Xq  will  be 
exceeded,  which  as  before  is  unsatisfactory.  The  cumulative  probability  density 
function  g(x)  =  G'(x)  is  related  to  the  initial  distribution  function  F(x),  and  the  modal  or 
most  probable  value  x^^  of  g(x)  is  in  fact  the  distribution  of  Pln  for  the  case  where 
Pli  =  l/N(Ochi,  1981),  as  shown  in  figure  7. 

If  the  initial  probability  density  fxmction  f(x)  is  assumed  to  be  Rayleigh  distributed, 
then  it  is  shown  that 


x„  =  ln(N)  (34) 

where  is  the  variance  associated  with  the  Rayleigh  distribution  f(x).  Ochi  (1981) 
suggests  the  use  of  a  risk  parameter  7  =  Pl^  to  be  applied  to  the  number  of  lifetime 
loading  cycles  N,  so  that  a  more  conservative  estimate  of  the  extreme  value  for  design 
purposes  x^  is  given  by 


Xq  =  V2ct"  ln(N/Y)  (35) 

which  combining  (34)  and  (35)  gives 
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If  for  example  7=0.01,  then  equation  provides  (36)  the  design  extreme  value  which  is 
not  expected  to  be  exceeded  with  99%  certainty.  It  can  be  seen  from  equation  (35)  that 
the  design  extreme  value  appears  to  be  more  sensitive  to  the  correct  evaluation  of 
than  it  is  to  either  N  or  7. 


Figure  7.  Long  Term  Probability  of  Exceedance  Based  Upon  Extreme  values 


5.2.4  Determination  of  the  Long  Term  Distribution  G(x) 

In  the  above  analysis  where  extrapolation  is  based  either  upon  the  modal  value or 
the  extreme  values,  it  is  necessary  to  determine  the  cumulative  long  term  probability 
distribution  G(x)  in  order  to  be  able  to  determine  the  design  load  for  the  ship 

The  wave  loading  of  a  ship  is  dependent  not  only  on  the  sea  conditions,  but  also  the 
ship’s  speed,  relative  heading  and  loading  condition,  and  throughout  the  life  of  the 
ship  there  will  be  a  certain  probability  of  occurrence  associated  with  any  particular 
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value,  or  range  of  values,  of  each  of  these  parameters.  The  short  term  probability 
calculated  from  trials  results  in  accordance  with  equation  (27)  relates  only  to  one 
particular  set  of  conditions,  which  will  occur  with  a  particular  probability  during  the 
life  of  the  ship.  Therefore  in  order  to  use  the  trials  information  to  determine  the  long 
term  distribution  it  must  be  weighted  according  to  these  probabilities.  Obviously  a 
large  amoimt  of  data  and  calculations  are  required  to  calculate  the  desired  probability 
and  to  simplify  the  problem  each  parameter  is  generally  considered  to  be  independent, 
although  this  is  not  strictly  correct.  The  long  term  cumulative  probability  may  then  be 
calculated  by  the  summation  of  the  short  term  probabilities  weighted  by  the  relative 
frequency  of  sea  state  and  the  relative  frequency  of  the  ship's  operating  state.  That  is 
the  long  term  probability  distribution  G(x)  is  given  by 


o2jC« 


G(x>Xo)  =  JJJJJ  Psi(x>Xo:V„X)/,Hv3.TJ.Pss(Hv3.TJ.Po(V„V,A).dV,.dH^3.dT,.(iX|/.dA 

(37) 


0  0  000 


where 

•  Pgi  is  the  short  term  probability  of  exceedance  and  is  dependent  upon  ship’s  speed, 
relative  heading  and  sea  state,  characterised  by  sigruficant  wave  height  and  zero 
crossing  period. 

•  Pgs  is  the  probability  of  encountering  a  particular  sea  state  characterised  by 
significant  wave  height  and  zero  crossing  period  and  which  is  determined  by 
voyage  analysis  and  reference  to  wave  atlases. 

•  Pq  is  the  probability  of  the  ship's  operating  condition  characterised  by  ship's  speed 
and  relative  wave  heading  (which  are  themselves  a  function  of  sea  state  although 
this  is  often  ignored  )  and  displacement.  (Note:  The  effect  of  ship's  speed  on 
bending  moment  is  often  quite  small  (Lewis  and  Zubaly,  1981)  and  is  often  ignored 
by  assuming  the  highest  possible  speed  for  a  given  set  of  conditions.) 

Although  equation  (37)  is  shown  as  a  continuous  integral,  in  practice  only  a  limited 
number  of  different  combinations  will  be  considered  as  the  level  of  calculation 
becomes  excessive.  The  probability  of  sea  state  conditions  may  be  determined  from 
wave  atlases,  which  use  observations  of  wave  height  and  period  (such  as  Hogben  et  al, 
1986)  to  predict  the  number  of  waves  per  1000  within  particular  wave  height  and 
period  bands  (see  section  2.5). 

The  probability  of  the  ship's  operating  condition  must  consider  the  distribution  of 
ship's  displacement,  ship's  speed  and  relative  direction  of  the  waves  compared  to  the 
ship's  heading.  This  is  determined  by  voyage  analysis  which  considers  the  time  that  a 
ship  will  spend  in  a  particular  area,  it's  speed,  heading,  displacement  and  the  relative 
seaway  direction.  For  specific  purpose  ships  such  as  bulk  carriers  the  response  to  a 
particular  seaway  will  vary  with  loading  condition,  and  so  for  such  ships  and  others 
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designed  for  a  particular  route  there  will  be  a  bias  towards  particular  speeds  and 
displacements  for  particular  legs  of  the  journey,  and  also  probably  towards  particular 
relative  seaway  directions.  For  more  general  purpose  and  Navy  ships  it  is  generally 
assiuned  that  each  relative  wave  direction  will  occur  with  equal  probability  so  this 
element  of  the  calculatioixs  may  be  omitted,  although  the  probability  of  speeds  and 
displacements  needs  to  be  considered. 


6.  Closing  Remarks 

The  study  of  waves  and  wave  induced  loading  is  a  complex  task  as  ocean  waves  are  a 
random  phenomenon  where  the  waves  form  an  irregular  pattern  of  wave  heights, 
lengths,  directions  and  frequencies.  Simplifying  techniques  such  as  spectral  and 
statistical  analyses  are  used  so  as  to  be  able  to  quantify  wave  properties  and  the  ship's 
responses  to  them.  Some  of  the  relevant  theories  have  been  presented  in  this  paper 
along  with  the  standardised  methods  with  which  they  may  be  applied  in  the  analysis 
of  ships’  structural  responses  to  wave  loading. 

The  relationship  between  actual  and  encoimter  wave  frequencies  is  discussed  and  the 
effect  of  relative  wave  direction  upon  this  relationship  is  highlighted.  The  discussion 
also  includes  the  development  and  use  of  wave  spectra,  the  different  forms  that  these 
spectra  may  take,  the  transformation  of  spectra  between  frequency  domains  and  the 
use  of  directional  spectra  to  account  for  wave  spreading.  The  typical  parameters  which 
may  be  used  to  quantify  waves  and  their  relationship  to  the  wave  spectra  are 
tabulated.  Idealised  wave  spectra  used  for  design  and  analysis  purposes  are  presented 
and  their  applicability  is  discussed. 

The  ship's  response  to  waves  is  often  assumed  to  be  linear  and  is  expressed  in  terms  of 
Response  Amplitude  Operators.  The  RAO  is  a  frequency  dependent  function  which 
provides  the  ratio  of  the  response  to  the  input  spectrum  over  the  frequency  range  of 
interest.  The  development  of  RAOs  requires  particular  attention  to  the  relationship 
between  encoimter  frequency  and  actual  wave  frequency  as  there  is  not  always  a  one- 
to-one  relationship  between  the  two. 

The  long  term  response  of  a  ship  to  wave  loading  over  a  wide  range  of  operating 
conditions  and  environments  will  ultimately  determine  the  fatigue  performance  of  the 
vessel.  As  the  sea  state  and  ship’s  operating  condition  are  continuously  changing  it  is 
necessary  to  make  long  term  fatigue  predictions  based  upon  short  term  observations. 
The  general  statistical  approach  to  short  term  and  long  term  use  of  such  data  is 
presented  and  the  applicability  of  particular  methods  is  discussed. 
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Nomenclature 


A 

constant  used  with  standard  wave  spectra 

B 

constant  used  with  standard  wave  spectra 

c 

wave  celerity 

m 

cosine  spreading  factor 

D 

variable  used  in  determination  of  spreading  factor  C 

e 

natural  logarithm,  2.71828 

E[] 

variance 

VEh 

modal  wave  height 

f 

cyclic  frequency 

f(x),  g(x) 

probability  density  function 

F(x),G(x) 

cumulative  probability  distribution 

g 

acceleration  due  to  gravity 

KVEh) 

long  term  probability  density  function  of  modal  wave  height  based  upon 
measured  values 

H 

wave  height 

H 

mean  wave  height 

H]/3 

significant  wave  height 

Hl/n 

average  height  of  1/n'th  highest  waves 

Hs 

visual  estimate  of  swell  wave  height 

Hv 

visual  estimate  of  wave  height 

Hw 

visual  estimate  of  wind  generated  wave  height 

Hgt?  (®e) 

transfer  fimction  of  response  with  wave  energy  spectrum 

m 

positive  integer 

mo 

zeroth  moment  of  the  spectral  area  about  the  origin 

mn 

n'th  moment  of  the  spectral  area  about  the  origin 

N 

number  of  cycles 

n 

positive  integer 
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n 


■Ll 


Po 

Psi 

Pss 

P(x,VEh) 

RAO 

S(co) 


S^((0,<t>) 

T 

T 

Tp 

Tp 

T. 

Tz 

Tz 

Vs 

Vw 

W(/) 


X 

Xd 

{Xt} 


Xi 


a 

tto 

P 

m) 


average  number  of  zero  upcrossings  per  unit  time 

average  lifetime  probability  of  exceedance  in  one  cycle  (modal  value 
formulation),  long  term  probability  of  exceedance  in  one  cycle  (extreme  value 
formulation) 

design  lifetime  probability  of  exceedance  (modal  value  formulation),long 
term  probability  of  exceedance  in  N  cycles  (extreme  value  formulation) 

probability  of  the  ship's  operating  condition 

short  term  probability  of  exceedance  in  one  cycle 

probability  of  encountering  a  particular  sea  state 

joint  long  term  probability  density  function  of  x  with  Ve  h 
Respoi:\se  Amplitude  Operator 

energy  spectrum  in  terms  of  circular  frequency  (subscript  denotes  measured 
parameter) 

directional  wave  energy  spectrum  m  terms  of  circular  frequency 
period  of  time 

"characteristic  period",  see  section  3.2 
wave  period  between  adjacent  peaks 
mean  wave  period 

visual  estimate  of  wave  period 

zero  crossing  period 
mean  zero  crossing  period 
velocity  of  the  ship 
wind  velocity 

energy  spectrum  in  terms  of  cyclic  frequency  (subscript  denotes  measured 
parameter) 

an  arbitrary  parameter 

estimate  of  the  extreme  value  for  design  purposes 
discrete  time  series 
extreme  response 

empirically  derived  coefficient  for  standard  wave  spectra 
wave  slope  amplitude 

empirically  derived  coefficient  for  standard  wave  spectra 
the  gamma  function 
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7 

M 

5(|) 

6(0 

z 

c 

Co 

Co 

c. 

X 

11 

C 

71 

P 

a 

<!> 

¥ 

CO 

CO 

COg 

®emax 

®k 


risk  parameter  for  long  term  probabilities 

discrete  time  interval,  inverse  of  digitising  frequency  / 

directional  spreading  band 

circular  frequency  interval 
bandwidth  parameter 
wave  elevation 
wave  amplitude 
mean  wave  amplitude 

most  probable  (modal)  wave  amplitude 

a  shape  parameter  used  with  the  Ochi  6-Parameter  Spectrum 

maximum(Xi,X2,...,Xn) 

the  mean  value  of  measured  a  parameter  x 

non-dimensional  parameter  x/  a 

pi,  3.14159 

density  of  sea  water 

the  standard  deviation  of  measured  a  parameter  about  it's  mean 
spreading  angle 

relative  heading  angle  between  ship's  heading  and  direction  of  waves 
circular  wave  frequency 
average  circular  frequency 

encounter  frequency  at  which  the  ship  encoxmters  waves 
maximum  encounter  frequency  for  given  conditions 
discrete  circular  wave  frequency 
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